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ABSTRACT 

Thermal emission from a rotating, supermassive star will cause the configuration to contract 
slowly and spin up. If internal viscosity and magnetic fields are sufficiently weak, the contract- 
ing star will rotate differentially. For each of the six initial angular momentum distributions 
considered, a cooling supermassive star will likely encounter the dynamical bar mode instability, 
which may trigger the growth of nonaxisymmetric bars. This scenario will lead to the generation 
of long wavelength gravitational waves, which could be detectable by future space-based laser 
interferometers like LISA. 

Subject headings: gravitation — instabilities — stars: evolution — stars: rotation 



1. Introduction 

The existence of supermassive black holes 
(SMBHs) and the theory of their role as fuel 
sources for active galactic nuclei and quasars are 
supported by a growing body of solid observa- 
tional evidence (see, e.g., the reviews of Rees 1998 
and Macchetto 1999, and references therein). For 
example, VLBI observations of the Keplerian mo- 
tion of masing knots orbiting the nucleus of NGC 
4258 indicate the presence of a central object of 
mass ~ 3.6 x 10 7 M© and radius < 13 pc (Miyoshi 
et al. 1995). The most reasonable explanation is 
that this extremely compact object is a black hole. 
High resolution HST observations of matter orbit- 
ing the nuclei of M87 and NGC 4151 also show 
nearly Keplerian motion around SMBHs (Mac- 
chetto et al. 1997; Winge et al. 1999). In addi- 
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tion, reverberation mapping and photoionization 
methods indicate the presence of SMBHs in the 
centers of 17 Seyfert quasars (Wandel, Peterson, 
and Malkan 1999). 

However, the nature of the progenitors of 
SMBHs is rather uncertain (see Rees 1984, 1999 
for overviews). SMBH formation scenarios that 
involve the stellar dynamics of dense clusters and 
the hydrodynamics of supermassive stars have 
both been proposed. The suggested stellar dy- 
namical routes focus on the evolution of dense 
star clusters. These clusters are formed via a con- 
glomeration of stars, produced by fragmentation 
of the primordial gas. In one such cluster scenario, 
massive stars form via stellar collisions and merg- 
ers in the cluster and then evolve into stellar-mass 
black holes. The merger of these holes, as they 
grow and settle to the cluster center, leads to the 
build up of one or more SMBHs (see, e.g., Quin- 
lan and Shapiro 1990, Lee 1995, and references 
therein). Alternatively, dense clusters composed 



of compact stars may become dynamically un- 
stable to relativistic core collapse and thus may 
give birth to SMBHs (Zel'dovich and Podurets 
1965; Shapiro and Tcukolsky 1985a,b; Quinlan 
and Shapiro 1987). 

The proposed hydrodynamical routes to SMBH 
formation center on the evolution of supermassive 
stars (SMSs). SMSs may contract directly out 
of the primordial gas or they may develop from 
a supermassive gas cloud that has been built up 
from the fragments of stellar collisions in clus- 
ters (Sanders 1970; Begclman and Rees 1978; 
Haehnelt, Natarajan, and Rees 1998; however, see 
also Loeb and Rasio 1994; Abel, Bryan, and Nor- 
man 2000; Bromm, Coppi, and Larson 1999). The 
evolution of SMSs will ultimately lead to relativis- 
tic dynamical collapse (Iben 1963; Chandrasekhar 
1964a, b) and thus possibly to the formation of 
SMBHs. 

In addition to their relevance to AGNs and 
quasars, SMBHs and their formation events are 
likely candidates for detection by proposed space- 
based gravitational wave detectors, like the Laser 
Interferometer Space Antenna (LISA) (Folkncr 
1998). LISA would be very sensitive to the long 
wavelength, low frequency gravitational radiation 
that supermassive objects are expected to emit 
(Thornc and Braginsky 1976; Thornc 1995). For 
example, LISA may be able to detect the collapse 
of a SMS to a SMBH or the coalescence of two 
SMBHs (LISA Prc-Phase A report 1995). The 
rates of such events are largely unknown due to 
the uncertainty in the formation mechanism of 
SMBHs. 

Identifying the scenario by which SMBHs form 
is hence of fundamental importance to a num- 
ber of areas of astrophysics. This work exam- 
ines one of the possible progenitors of SMBHs dis- 
cussed above, SMSs. The structure and evolution 
of SMSs have been explored by many authors (see, 
e.g., Zel'dovich and Novikov 1971 and Shapiro and 
Tcukolsky 1983, and references therein). In this 
paper, we focus specifically on the quasistatic evo- 
lution of SMSs as they contract and cool due to 
thermal emission. 

Two related factors that affect the cooling 
phase of a SMS are the strengths of its viscos- 
ity and magnetic fields and the nature of its an- 
gular momentum distribution. If strong enough, 
viscosity will enforce uniform rotation throughout 



the star as it contracts. However, the strength of 
molecular viscosity is unlikely to be large enough 
to accomplish this, so turbulent viscosity or mag- 
netic fields must provide the necessary dissipation 
(Zel'dovich and Novikov 1971; Spruit 1999a,b; 
Baumgarte and Shapiro 1999). However, it is 
not known whether these agents are sufficiently 
strong to maintain uniform rotation in SMSs as 
they evolve. 

If a SMS is rotating prior to its contraction 
phase, conservation of angular momentum re- 
quires the star to spin up and therefore increase 
its ratio of kinetic energy T to gravitational po- 
tential energy W, (3 — T/\W\. If viscosity or 
magnetic fields are able to maintain uniform rota- 
tion, the star will be spun up to its mass-shedding 
limit f3 s hed — 8.99 x 10~ 3 . The value is small be- 
cause a SMS is dominated by radiation pressure 
and thus is well represented by an n=3 polytropic 
equation of state (see §2.1.1). As a result, it is 
quite centrally and cannot support much rotation 
when spinning uniformly before mass at the equa- 
tor is driven into Keplcrian orbit. Thereafter, it 
will evolve along a mass-shedding sequence, los- 
ing both mass and angular momentum, and will 
contract to the point of relativistic instability and 
collapse. Baumgarte and Shapiro (1999) have 
recently investigated this evolutionary scenario 
in detail and have described how it may lead to 
SMBH formation. Specifically, they identified the 
critical configuration at the onset of the relativis- 
tic radial instability. They also demonstrated that 
the nondimcnsional (in units where G=c=l) ra- 
tios T/|W|; R/M, where R is the radius and M 
is the mass; and J/M 2 , where J is the angular 
momentum, are universal numbers (i.e., they do 
not depend on M, R, or J for the crititcal config- 
uration or on the history of the star). 

This work investigates the extreme opposite 
regime of SMS cooling evolution, the low- viscosity 
and low-magnetic field limits. In this case, as the 
star contracts and spins up, neither viscosity nor 
magnetic fields will be able to enforce uniform ro- 
tation (even if the star initially rotates uniformly) . 
The outcome of the evolution will then depend on 
the star's initial angular momentum distribution 
(or rotation law) , as the star will conserve its angu- 
lar momentum on cylindrical shells as it contracts. 
One possible outcome is that the star will again be 
spun-up to the mass-shedding limit (at a value of 
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fished higher than that for uniform rotation) . The 
subsequent evolution is not clear. It is likely that 
mass and angular momentum will be ejected from 
the equator as the star undergoes further cooling. 
It is possible that it will continue to evolve in qua- 
sistationary equilibrium until reaching the onset 
of dynamical collapse. If it implodes to a SMBH, 
the matter is likely to emit an appreciable burst 
of gravitational radiation just prior to entering the 
event horizon. 

An alternative outcome of low-viscosity, low- 
magnetic field contraction of a differentially ro- 
tating SMS is that it will encounter the dynami- 
cal "bar mode" instability at a critical value (ib ar 
prior to reaching fished- Such global rotational 
instabilities in fluids arise from nonaxisymmetric 
modes, which vary like e ±M ™^, where m=2 is the 
bar mode. As the bar mode develops, angular mo- 
mentum and mass will be transported outwards 
(Tohline, Duriscn, and McCollough 1985; Pickett, 
Durisen, and Davis 1996; Imamura, Durisen, and 
Pickett 2000; New, Centrella, and Tohline 2000; 
Brown 2000). 

If enough angular momentum is removed, col- 
lapse to a SMBH could ensue. A SMS undergoing 
a bar mode distortion would be a strong source 
of quasi-periodic gravitational radiation because 
of its highly nonaxisymmetric structure. We note 
that a bar might form during catastrophic collapse 
in any scenario, if /3(oc 1/R) grows above the crit- 
ical value for the dynamical bar mode instability. 

The outcome of the cooling and the resul- 
tant spin-up of a low-viscosity SMS depends crit- 
ically on its initial angular momentum distribu- 
tion. However, the nature of the angular momen- 
tum distribution of SMSs is largely unknown. The 
purpose of this paper is to examine the evolution 
of SMS models by assuming several different ini- 
tial rotation laws to determine, for example, if any 
of them might cause the star to encounter the dy- 
namical bar mode, prior to mass-shedding and col- 
lapse (i.e., to determine if (ibar < fished)- 

To make this determination, we analyze several 
Newtonian equilibrium sequences of SMS models. 
The models along each sequence are constrained 
to have the identical mass M, angular momentum 
J, and angular momentum distribution but de- 
creasing entropy, due to cooling. Such a sequence 
mimicks the quasistatic evolutionary track of a 
SMS as it slowly cools and contracts. We find that 



along the sequence, the stars flatten as measured 
by the ratio of polar radius R p to equatorial ra- 
dius R e . We examine each sequence to determine 
whether the limits fished and/or fiber are reached. 
Note that 3D hydrodynamical simulations are nec- 
essary to determine precisely fibar for each se- 
quence. However, previous linear and nonlinear 
stability analyses indicate that 0.14 < fib ar ^ 
0.27 (Chandrasekhar 1969; Shapiro and Teukolsky 
1983; Durisen and Tohline 1985; Imamura, Fried- 
man, and Durisen 1985; Managan 1985; Hachisu, 
Tohline, and Eriguchi 1988; Tohline and Hachisu 
1990; Pickett, Durisen, and Davis 1996; Toman 
ct al. 1998; Imamura, Durisen, and Pickett 2000). 
See §4 for further discussion of these stability lim- 
its. 

We note that secular version of the bar insta- 
bility, which sets in at a lower value of fi, is not 
relevant for the scenario described here. The rea- 
son for this is that the viscosity that could drive 
the mode is assumed to be vanishingly small and, 
hence, the secular timescale is too long (Baum- 
garte and Shapiro 1999). The timescale for grav- 
itational radiation to drive the secular bar mode 
is also much longer than the evolution (cooling) 
timescale (Chandrasekhar 1970). 

In the low-viscosity and low-magnetic field lim- 
its, the angular momentum of mass shells on con- 
centric cylinders in a contracting star will be pre- 
served (Bodcnhcimer and Ostriker 1973; Tassoul 
1978). Thus we confine our discussion to equi- 
librium sequences constructed with rotation laws 
that conserve the angular momentum distribution 
in this manner. The sequences examined here in- 
clude one we construct below with the n'=3 rota- 
tion law (see §2.2). This sequence represents the 
cooling evolution of a SMS that is nearly spher- 
ical and rotating uniformly at formation prior to 
its contraction phase. It is the low-viscosity, low- 
field analog of the uniformly rotating sequence 
considered in Baumgartc & Shapiro (1999); in the 
present case, differential rotation will ensue once 
cooling gets underway. Because the initial rota- 
tion profiles of SMSs are unknown, we also exam- 
ine alternative sequences constructed by Hachisu, 
Tohline, and Eriguchi (1988; hereafter HTE) that 
describe the evolution of SMSs with nonuniform 
initial rotation laws. 

The remainder of this paper is organized as 
follows. In §2, we describe the numerical meth- 
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ods used to construct our n'=3 equilibrium se- 
quence and the sequences constructed by HTE. 
We present our n'=3 sequence and summarize the 
sequences constructed by HTE in §3. In §4, we 
outline possible evolutionary scenarios for SMSs, 
based on the sequences described in §3. We discuss 
these results in §5. 

2. Numerical Methods for SMS Model 
Construction 

We have constructed an equilibrium sequence of 
SMS models, with a rotation law appropriate for 
SMSs that rotate uniformly prior to contraction. 
The individual models along the sequence have the 
same rotation law, M, and J, but decreasing en- 
tropy and axis ratios R p /R e . Each sequence is 
a quasistatic approximation to the evolution of 
a SMS that is contracting due to thermal emis- 
sion. Modeling this phase of the evolution of SMSs 
with such equilibrium sequences is appropriate as 
the cooling timescale is longer than the hydrody- 
namic timescale (for M < 10 13 M Q , Baumgarte 
and Shapiro 1999, and references therein) 

2.1. The Self-Consistent Field Method 

The structure of a fluid rotating about the z 
axis with constant angular velocity f2 = Q(zu), 
where w is the distance from the rotation axis, is 
described by the following expression: 

-WP + V$ + hlVV(m) = 0, (1) 
P 

where P is the pressure, $ is the gravitational po- 
tential, ^(w) — — (1//iq) J fl 2 (w)zu dm is the cen- 
trifugal potential, and ho is a constant. Such a 
fluid is said to be in hydrostatic equilibrium be- 
cause the forces due to its pressure and to its grav- 
itational and centrifugal potentials are in balance. 

The method we have used to construct SMSs in 
hydrostatic equilibrium is Hachisu's grid based, it- 
erative, axisymmetric Self-Consistent Field (HSCF) 
technique (Hachisu 1986, see also New 1996). 

2.1.1. Equation of state 

The HSCF method requires the fluid to have 
a barotropic equation of state (EOS) P — P(p). 
All of the SMS models we have constructed have 
a polytropic EOS for which 

P = Kp 1+1 ' n , (2) 



where n is the polytropic index and K is the en- 
tropy constant. The enthalpy of a polytrope is 



H = (n + l)Kp 



l/r 



(3) 



The structure of a SMS is well represented by an 
n=3 polytrope (see §17.2 of Shapiro and Teukol- 
sky 1983). This is because SMSs are dominated 
by thermal radiation pressure (for M > 10 6 M Q , 
Zcl'dovich and Novikov 1971; Fuller, Woosley, and 
Weaver 1986) and are convective with constant 
entropy per baryon throughout (see Loeb and 
Rasio 1994 for a simple proof of this property). 
The value of K is determined by the entropy per 
baryon s in the star and is given by 



K 



a ( 3s \ 
3 \4m H a J 



4/3 



(4) 



where a is the radiation density constant and tuh 
is the mass of a hydrogen atom (see eq. 17.2.6 in 
Shapiro and Teukolsky 1983). The constant s can 
be expressed in terms of the matter temperature 
T m , baryon denisty rib, an d a: 



3 n b 



(5) 



For polytropic EOSs, the sequences constructed 
with Hachisu's SCF method are parameterized 
solely by the axis ratio R p /R e . The method re- 
quires that the grid cell locations of the two sur- 
face points, A = R e and B = R p , be specified. 

2.1.2. Rotation Laws 

In the low- viscosity limit, the angular momen- 
tum of a contracting star will be preserved on 
cylindrical mass shells (Bodenhcimcr and Ostriker 
1973; Tassoul 1978). Thus the rotation laws used 
in sequences representative of such evolution must 
enforce this conservation. The n'=3 law and each 
of the laws employed in the sequences constructed 
by HTE are defined such that the specific angu- 
lar momentum profile j(m) is identical for each 
model constructed. Here m is the mass interior 
to cylindrical radius w. Because j is the angu- 
lar momentum per unit mass and is constant on 
cylindrical surfaces, a j{m) law does ensure that 
the angular momentum distribution on cylindrical 
mass shells will be preserved from model to model 
along the sequence. 
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2.1.3. Converged Configurations 

The quality of the converged configurations can 
be estimated in terms of the virial error VE. The 
VE provides a measure of how well the energy is 
balanced and is defined as 



VE = 2T + W + 3 / PdV. 



I 



(6) 



where V is the volume of the model. The quantity 
VE is zero in strict equilibrium. The VEs for 
models on this sequence range from ~ 10 -2 (for 
the very smallest axis ratios) to ~ 1CP 6 , where the 
individual terms in equation (9) are of order unity 
(see Table 1). 

Our axisymmetric equilibrium configurations 
were constructed on uniform, 1024 x 1024 (zu, z) 
grids. Equatorial symmetry through the z~0 
plane was assumed; this symmetry condition is 
valid for systems with and barotropic EOSs 

(Tassoul 1978). 

Note that initially we had difficulty obtaining 
converged models on the n=n'=3 sequence (see 
§2.2, 3.1) with axis ratios < 0.400. For these mod- 
els, the HSCF iterations oscillated between two or 
more states without converging. In order to obtain 
converged models with these axis ratios, we im- 
plemented a modification to the HSCF algorithm 
suggested in Pickett, Durisen, and Davis (1996). 
Instead of using from iteration i as input for the 
next iteration step i + 1, a combination of pi and 
the old density Pi-i is used ("undcr-relaxation"): 



Pi = (! - &)Pi + $Pi 



(7) 



where S is a constant (in this work, 0.75 < S < 
0.99). 

2.2. Scale of Models 

As mentioned above, polytropic sequences con- 
structed with the HSCF code are parameterized 
soley by the axis ratio R p /R e . That is, for a given 
EOS and rotation law, the choice of R p /R e pro- 
duces a unique model. A dimensional scale (i.e., 
mass, radius, etc.) for this model may be cho- 
sen after its construction. The individual models 
along a sequence representative of the cooling evo- 
lution of a SMS have the same M and J. As will 
be demonstrated below, the choice of M and J 
for a sequence sets the dimensional scale for its 
individual models. 



The HSCF computations are carried out in di- 
mensionless form, with all quantities normalized 
in terms of G, R e , and p max . Here p ma x is the 
maximum density of the model. In HSCF units, 
M and J are normalized according to: 

M 

M = (8a) 

1 t e f J max 

J = i J 2 ■ (8b) 
GiR\p?nax 

Here, and in what follows, carets denote normal- 
ized quantities. 

The HSCF method produces unique values of 
the normalized quantites (e.g., M, J, etc.) for 
each input axis ratio R p /R e . After a converged 
model has been constructed, the normalized quan- 
tities are known. If physical values of M and J 
are then chosen to specify the sequence, equations 
(8) can be used to compute R e and p max for each 
model along the sequence: 



R< 




Pmax — 



MRl 



(9) 



(10) 



Once R e and p max are computed, they can be 
used to convert normalized quantities to phys- 
ical units. For example, the definition of H 
from equation (3) and its normalization, H = 
H/(GRlp m ax), can be used to determine a physi- 
cal value for K: 

1 



K 



1 + n 



GHrnaxRlp^aJ 71 . (11) 



Note that when n=3, equation (10) leads to the 
following simplification of equation (11): 



4 max \ M 



(12) 



This is the well known result that the scale of K 
depends solely on the choice of M, for n='i poly- 
tropes. 

3. Equilibrium Sequences 

In this section, we summarize the properties of 
n=3 polytropic equilibrium sequences. The n'=3 
sequence constructed by the present authors is dis- 
cussed in §3.1 and the sequences constructed by 
HTE are discussed in §3.2. 
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3.1. n'=3 

As mentioned in §1, the n'=3 law corresponds 
to the rotation profile of a uniformly rotating, 
spherical n—3 polytrope. In the low-viscosity 
limit, a sequence constructed with this law is rep- 
resentative of the evolution of a SMS that rotates 
uniformly prior to its cooling/contraction phase. 
Thus this sequence begins with the same starting 
model as that used in Baumgarte and Shapiro's 
study of SMS evolution in the high-viscosity limit 
(Baumgarte and Shapiro 1999). The n'=3 specific 
angular momentum distribution is given by (Bo- 
denheimer and Ostrikcr 1973): 

.. . / m{w)\ a2 / m (vj)\ a3 
j(m) = a 1 +a 2 (l-^ r >) +^^1 , 

(13) 

where M is the total mass of the system, and the 
numerically determined constants are a\ = 13.27, 
a 2 = 163.3, a 3 = —176.5, a 2 = 0.2353, and a 3 = 
0.2222. 

In Table 1, we list the properties of selected con- 
figurations constructed for the n'=3 rotation law. 
The data given (in HSCF normalized units, see 
§2.2) include the axis ratio R p /R e , polytropic en- 
tropy constant K (normalized to the entropy con- 
stant of the initial spherical model Kq) 1 mass M, 
angular momentum J, maximum enthalpy H m ax, 
total energy E = E/GR e 5 p max 2 , ratio of kinetic 
to potential energy /3, ratio of the centrifugal force 
F c to the gravitational force F g on the equator, 
and the virial error VE. When the quantity 
(F c /F g ) eqi equals unity, the mass shedding limit 
fished has been reached. 

Figure 1 displays R p /R e , [3, (F c /F g ) eq , and E 
as a function of Kj Kq for the individual models on 
this sequence. Note that K, which is a measure of 
the specific entropy, decreases monotonically with 
Rp/Re along the sequence. (The ratio K/K can 
be related to time once the cooling law is specified; 
this has been treated by Baumgarte and Shapiro 
1999 for the case of uniform rotation.) The energy 
E also decreases monotonically. Density contour 
plots of selected sequence models are shown in Fig- 
ure 2. 

This sequence terminates due to mass shed- 
ding. The last spheroidal model we were able to 
construct had an axis ratio R p /R e = 0.002 and 
(F c /F g )eq— 0-99. We were unable to construct a 
converged model with i? p /i? e =0.001, as the it- 



eration process oscillated between models with 
(F c /Fg) less than and greater than one. Thus 
we conclude that this sequence terminates due to 
mass shedding at or very near R p /R e ~ 0.001, with 
fished ~ 0.30. Note that there is no physical effect 
that forces the mass shedding limit to occur at an 
axis ratio of 0. In fact, these two points do not co- 
incide on equilibrium sequences constructed with 
other equations of state and/or rotation laws (see, 
e.g, Hachisu 1986). 

Note that no toroidal branch of this sequence 
exists. The HSCF code can construct tori if the in- 
ner radius of the torus Ri is specified instead of R. p . 
All of the torioidal configurations we attempted to 
construct did not converge for this sequence (and 
had {F c /Fg) eq > 1 during the iteration process). 

Recall that the dynamical bar mode instabil- 
ity is likely to set in at f3b ar < 0.27. The model 
with (3 ~ /?b ar =0.27 on the n'=3 sequence has 
i?p/i? e =0.004 and K/K Q =0.S81. 

In order to estimate the possible physical 
scales of SMS models on this sequence, we 
must choose the constant values M and J (see 
§2.2). For example, say the initial nearly spher- 
ical star has M=10 6 M Q , (i? e ) =10 17 cm, and 
Po = (T/|VF|)o=10- 5 , where a subscript de- 
notes the value is for the (nearly) spherical star. 
The value of (R e )o is chosen such that the com- 
paction, (G/c 2 )M/ R, for the model that reaches 
~ f3bar—0-^7 is of the same order of magnitude 
as the uniformly rotating critical configuration of 
Baumgarte and Shapiro (1999) for gravitational 
collapse due to relativistic gravitation. For all 
larger values of the product f3o(R e )o , an evolved 
configuration with M=10 6 Mq will experience a 
dynamical bar instability prior to reaching the 
critical compaction for radial collapse. The value 
of /3o is quite arbitrary except for the slow rotation 
constraint, /?o "C 1, necessary to ensure that the 
initial cloud is spherical. Cosmological considera- 
tions and calculations of early structure formation 
will ultimately determine the likely distribution of 
parameters for primordial SMSs. 

Having chosen M , (R e )o and Po, one can exhibit 
the scaling of all of the remaining quantities de- 
rived below in terms of these choices, i.e. in terms 
of (M 6 ee M/1O 6 M ), (R 17 ee (i? e ) /10 17 cm), 
and (/?_ 8 ee/VIO- 5 ). 

The choice of M, (R e )o an d /3o determines 
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J. To see this, recall that for a uniformly ro- 
tating star J = (2T/) 1 / 2 , where I is the moment 
of inertia. On this sequence, |M / |o=8.96 x 10~ 3 
and 7 =5.82 x 1(T 3 . Thus, J = 3.23 x Kr 5 /? V 5 2 . 
Equations (8b) and (10) can then be used to con- 
vert J to J=1.09xl0 61 [/3_ 5 i2i 7 M|] 1 /2g Cm 2 s -i. 

These choices of M and J for the sequence can 
be used to dimensionalizc its individual models 
according to equations (9-12). For instance, the 
nearly spherical initial model has (p m ax)o= 2.57 x 
10~ 11 M 6 i? 17 3 gcm- 3 andKo=3.84xl0 18 Mg /3 g- 1 / 3 cm 3 
The model with (3 ~ /3b ar — 0.27 and axis ratio 
0.004 has # e =4.14 x 10 13 /3_ 5 i?i 7 cm, p max =1.27 x 
10 5 M 6 [/?_ 5 i?i 7 ]- 3 gcm- 3 , andK=1.46xl0 18 M 6 2/3 g- 1 / 3 
The model with axis ratio 0.002 has i? e =4.04 x 
10 13 /3_ 5 i?i 7 cm, p mQ:c =4.18xl0 5 M 6 [/3_5/i 17 ]- 3 gcm- 3 , 
and #=1.20 x 10 18 M 6 2/3 g- 1/3 cm 3 s" 2 . 

3.2. Hachisu, Tohline, and Eriguchi 1988 

Hachisu, Tohline, and Eriguchi (1988), here- 
after HTE, used the SCF method to construct 
n=3 polytropic sequences with the following pa- 
rameterized angular momentum distribution: 

j(m) = (l + q)(J/M)[l-(l-m(w)/M) 1/q ]. (14) 

Here, the index q specifies the rotation law. Note 
that the limiting case of q=0 corresponds to the 
j-constant rotation law. 

The j-constant (q=0) sequence HTE con- 
structed is entirely toroidal. That is, no spheroidal 
models exist with this rotation law (for any of the 
polytropic indices < n < 3 they investigated). 
Thus the j-constant sequence is not relevant to the 
discussion of the cooling phase of SMS evolution, 
which begins with a spheroidal configuration. 

HTE also constructed sequences with q=0.5, 
0.7, 1.0, 1.3, and 1.5. All of these sequences 
make continuous transitions between spheroidal 
and toroidal configurations at relatively high val- 
ues of Ptrans > 0.33, without encountering mass 
shedding. Note that if M and J are held constant, 
the entropy K and energy E decrease monotoni- 
cally along the spheriodal branch of each of these 
sequences. This confirms that they are suitable 
representations of the cooling evolution of SMSs. 

From HTE's Figure 7f, one can estimate that 
Kbar/K a ~ 0.4 for each of these sequences (where 
K bar is the value of K for a model with /3=0.27). 
Thus the onset of the bar mode instability for 



HTE's sequences occurs at approximately the 
same value of K/K as on the n'=3 sequence. 

4. Evolutionary Scenarios 

In the low- viscosity limit, the outcome of the 
cooling evolution and resultant spin up of a SMS 
depends critically on its initial angular momen- 
tum distribution. Because the nature of the angu- 
lar momentum distribution in SMSs is largely un- 
known, we have presented equilibrium sequences 
s of SMS models with several different rotation laws 
in §3. Here we outline possible scenarios for the 
outcome of the cooling evolution of SMSs with low 
cn ¥if?cosity, based on these sequences. 

As mentioned in §1, one possible outcome is 
that a SMS will spin up to a critical value f3b ar , 
at which it will encounter the dynamical bar mode 
instability. As the bar mode develops, angular mo- 
mentum and mass will be transported outwards. 
The angular momentum removed will hasten, to 
some degree, the eventual collapse of the SMS (and 
the subsequent possible formation of a SMBH). 
A SMS undergoing the bar mode is potentially a 
strong source of gravitational radiation because of 
its highly nonaxisymmetric structure and hence 
might be detectable by LISA. 

The nonlinear nature of the bar mode insta- 
bility necessitates the use of 3D hydrodynami- 
cal simulations to precisely determine fi^ar for 
each equilibrium sequence. However, previous lin- 
ear and nonlinear stability analyses indicate that 
0.14 < flbar ^ 0.27 for a wide range of polytropic 
indices and rotation laws (Shapiro and Teukolsky 
1983; Durisen and Tohline 1985; Imamura, Fried- 
man, and Durisen 1985; Managan 1985; Hachisu, 
Tohline, and Eriguchi 1988; Tohline and Hachisu 
1990; Pickett, Durisen, and Davis 1996; Toman 
et al. 1998; Imamura, Durisen, and Pickett 2000). 
The lower limit, 0.14, may be most appropriate 
for spheroids with off-center density maxima and 
tori (Hachisu, Tohline, and Eriguchi 1988; Tohline 
and Hachisu 1990). 

The nonlinear hydrodynamics study of Pickett, 
Durisen, and Davis (1996) indicates that 0bar may 
be less than 0.27 for rotation laws that place a 
significant amount of angular momentum in equa- 
torial mass elements. Their results predict that 
the m=2 stability limit may be less than ~ 0.20 
for models with the n'—3 rotation law. Their sim- 
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ulations, of n=1.5 polytropes, also suggest that 
a one-armed spiral, m—1 mode may become in- 
creasingly dominant over the to=2 mode as the 
equatorial concentration of angular momentum in- 
creases. Note that the grid resolutions used in 
Pickett, Duriscn, and Davis (1996) were likely not 
sufficient to accurately model the development of 
instabilities in models with these extreme differ- 
ential rotation laws (Toman et al. f998). How- 
ever, the results of Pickett, Durisen, and Davis 
(1996) and the linear and nonlinear stability anal- 
yses of Toman et al. (1998) confirm that 0.27 is 
likely to be an approximate upper limit to fibar, for 
a variety of polytropic indices and rotation laws. 
The analysis presented in the remainder of this 
manuscript assumes that the m=2 bar mode is 
the dominant mode and that fibar ^ 0.27. 

The new n'=3 sequence constructed by the 
present authors and the sequences constructed by 
HTE that are considered in this work indicate 
that a bar mode phase is probable for SMSs with 
a wide range of differential rotation laws. In- 
deed, HTE's q=0.1-1.5 (see equation 14) sequences 
do not terminate due to mass shedding. Each 
of these sequences makes a continuous transition 
from spheroidal to toroidal configurations at val- 
ues of Ptrans > 0.33 > 0.27 > (i bar . Thus, n=3 
models with these g-indexed laws would likely en- 
counter the bar mode as spheroids. 

The spheroidal n'—3 sequence terminates due 
to mass shedding at fished > 0.30. Because 
fished > 0.27, we expect that models with this ro- 
tation law become dynamically unstable to the bar 
mode near R p /R e ~ 0.004 (this model is marked 
with a solid dot in Figure 1), before reaching the 
mass shedding limit. 

Note that even if the actual value of (ibar is 
less than 0.27, the qualitative nature of our re- 
sults would not change. That is, the sequences 
examined here would still have models that are 
unstable to the bar mode. In addition, the quan- 
titative estimates of the characteristics of the grav- 
itational radiation emission presented in §5 would 
only change by a numerical factor of order 1-10, 
even if (ib ar were as low as 0.22. 

3D hydrodynamical simulations are needed to 
properly study the evolution of rapidly rotating 
stars that encounter the bar mode instability. To 
date, we are unaware of any published hydrody- 
namical studies of the bar mode in n=3.0 poly- 



tropes. However, many authors have simulated 
this instability in polytropes with lower n and var- 
ious rotation laws (e.g., Tohline, Durisen, and Mc- 
Collough 1985; Pickett, Durisen, and Davis 1996; 
Imamura, Duriscn, and Pickett 2000; New, Cen- 
trella, and Tohline 2000; Brown 2000). All of these 
simulations agree on the initial phase of the in- 
stability, during which the initially axisymmctric 
object deforms into a bar shape through which 
spiral arms are shed. Recent Newtonian simu- 
lations indicate that the long-term outcome of 
the bar instability in an n=1.5 polytrope (with 
an initial Maclaurin rotation law) is a persistent 
bar-like structure that emits gravitational radia- 
tion over many rotation periods (New, Centrella, 
and Tohline 2000; Brown 2000) (the persistent na- 
ture of the bar structure was recently confirmed 
with the post-Newtonian simulations of Saijo et 
al. 2000). Simulations relevant to the dynamical 
bar mode instability have recently been performed 
for n—1 polytropes in full 3D general relativity by 
Shibata, Baumgarte, and Shapiro (2000). 

5. Conclusions and Future Work 

The thermal emission of a rotating supermas- 
sive star causes it to contract and spin up. If the 
viscosity and internal magnetic fields are weak, the 
SMS will rotate differentially during its cooling 
and contraction phase. We have examined mod- 
els of the quasistatic cooling evolution of SMSs, 
represented by equilibrium sequences with six dif- 
ferential rotation laws. 

The evolution of a cooling SMS that is initially 
nearly spherical, with very low, uniform rotation 
can be represented by the n'=3 sequence we con- 
structed. This sequence terminates due to mass- 
shedding at fished > 0.30. This fished exceeds the 
value 0.27, which is a likely upper limit for the 
onset of the rotationally induced, dynamical bar 
instability. Thus it appears likely that a star with 
this rotation law would become unstable to the bar 
mode prior to reaching the mass shedding limit. 

We have also considered sequences constructed 
by HTE, in order to follow the evolution of SMSs 
with five different initial differential rotation laws. 
No mass-shedding limits exist on these sequences. 
Thus, a cooling SMS characterized by one of these 
rotation profiles will also encounter the dynamical 
bar mode. 
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This mode will transport mass and angular mo- 
mentum outward and thus may hasten the onset 
of collapse and possible SMBH formation (Tohlinc, 
Durisen, and McCollough 1985; Pickett, Durisen, 
and Davis 1996; Imamura, Durisen, and Pickett 
2000; New, Centrella, and Tohlinc 2000; Brown 
2000). Because of their nonaxisymmctric struc- 
ture, SMSs undergoing the bar mode could be 
strong sources of quasi-periodic gravitational radi- 
ation for detectors like LISA, even before collapse. 

The frequency of these quasiperiodic gravita- 
tional waves can be estimated from the expected 
bar rotation rate flbar- The model on our n'=3 se- 
quence with (3 ~ /?6 ar =0.27 has a central rotation 
rate fi c =1.02 x lO^Hz [M 6 /rfi?7/ 7 3 ] 1/2 . In previ- 
ous hydrodynamics simulations of the bar mode 
instability (New, Centrella, and Tohline 2000), 
ilbar was ~ 0.4O c . With this relation between 
ilbar and fl c , the gravitational wave frequency 
few can be estimated to be 



few = Zfbar = 2 



a 



bar 



0.4, 



2tt 



IT 

3 D-3l 



~ 1 x IQ- 2 Hz[M 6 /3IIR^}2. (15) 

For a SMS of 1O 6 M0, which begins as a slowly 
rotating star of radius 10 17 cm with /3o=10 -5 , this 
yields a frequency of 1 x 10 _2 Hz. This frequency 
is in the range in which LISA is expected to be 
most sensitive, 10~ 4 -lHz (see, e.g., Thorne 1995, 
Folkner 1998). The choice /3 =2 x \Q- i /R 17 is the 
maximum value of (3q for which /gvk (=10~ 4 Hz) 
is still in LISA's range of sensitivity. 2 

The strength of the gravitational wave signal 
can be estimated roughly to be 



h 



GQ 
c 4 d 



G MR barfbar 



4 x 10~ 15 



lGpc 



MipllRrf, (16) 



where Q is the second time derivative of the star's 
quadrupole moment and d is the distance, which 
we scale to 1 Gpc (the Hubble distance is ~ 3 
Gpc). Here we have used Rb ar ~ i? e =4.14 x 



2 We note that Peebles (1969) has discussed the spin-up of 
primordial masses by tidal torquing and finds that J ~ 
MB? . Naively scaling his results (derived for M ~ 10 11 Afg 
configurations) to our regime, we find that /3o ~ 10~ 4 , 
so that the resulting frequency of gravitational radiation 
would lie near the lower limit of LISA's range. 



10 13 /3_5i?i7 cm for the n'=3 model that reaches 
the point /3 ~ Pbar- 

As mentioned previously, recent simulations of 
the bar instability indicate that the outcome is 
a long-lived bar-like structure (New, Centrella, 
and Tohlinc 2000; Brown 2000; Saijo et al. 2000). 
The bar will decay on a secular timescale due 
to dissipative effects. For differentially rotating 
SMSs, the largest source of dissipation will be 
gravitational radiation. The gravitational radia- 
tion damping timescale tqw is approximately 

T 

T GW ^ JJe\ ' 

V dt I GW 

where T is the rotational kinetic energy and 
(dE I 'dt)aw is the rate at which gravitational radi- 
ation carries energy away from the system. Recall 
that T = The gravitational potential en- 

ergy \W\ ~ GM 2 /R bar . Thus, 



GM 2 
R 



(18) 



The radiation rate can be estimated as (Shapiro 
and Teukolsky 1983) 



dE 



G M 



dt lGW~ c 5 R 2 V 



(19) 



In this case the characteristic velocity of the sys- 
tem is v — {(3GM/R) 1 / 2 . Substitution of equations 
(18) and (19) into equation (17) yields 



R 



bar 



G 3 /3 2 M 3 

1 x 10VsM 6 - 3 [/3_5i?i7] 4 . (20) 



The number of cycles J\f for which the signal 
will persist is 

A/" = tgw/gw 

= 4 x 10 9 [M 6 /3_ 5 i?i7r 5/2 . (21) 

The quasiperiodicity of such a signal will assist in 
its detection (Schutz 1997). 

The fraction of the mass (dM /M)gw radiated 
via gravitational radiation over the interval tqw 
can be estimated as 



dM \ 
~M)gw 



,dE, 



1 



TGW\ 



dt J GW Mc 2 



G f3M 
^~R 

1 x 10- 3 M 6 [/3_ 5 i? 17 ]- 1 .(22) 
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It appears that SMSs with a wide range of dif- 
ferential rotation laws may become unstable to the 
development of the dynamical bar mode. How- 
ever, linear and nonlinear stability analyses are 
needed to precisely determine the onset of the 
m=2 dynamical instability for the equilibrium se- 
quences considered here. Such analyses are also 
needed to determine the relative importance of 
various unstable nonaxisymmetric modes in SMS 
models (as a dominant m—1 mode would change 
the characteristics of the gravitational radiation 
emission). In addition, three-dimensional hydro- 
clynamical simulations are necessary to study the 
mass and angular momentum redistribution in- 
duced by these nonaxisymmetric instabilities, to 
compute the gravitational waveforms emitted, and 
to determine the final fate of the star. 

Future investigations involving hydrodynamical 
simulations of SMS models would benefit from im- 
proved knowledge of initial conditions, such as an 
appropriate value for (3 - These appropriate ini- 
tial conditions could be determined from studies 
of large-scale structure and cosmology. 

Most importantly, studies are required to as- 
sess the roles of viscosity and magnetic fields in 
rotating SMSs to judge whether they are sufficient 
to drive these configurations to uniform rotation 
prior to bar instability. 

We thank Joel Tohline for helpful conversa- 
tions. This work has been supported in part by 
NSF Grants AST 96-18524 and PHY 99-02833 and 
NASA Grants NAG5-7152 and NAG5-8418 to the 
University of Illinois at Urbana- Champaign. A 
portion of this work was performed under aus- 
pices of the U.S. Department of Energy by Los 
Alamos National Laboratory under contract W- 
7405-ENG-36. 
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Fig. 1. — Properties of individual models on the 
n'=3 equilibrium sequence are displayed (IS (V func- 
tion of the ratio K/Ko, where Kq is the en- 
tropy constant for the initial, spherical model: (a) 
Rp/ R e , axis ratio, (b) /3, ratio of kinetic to grav- 
itational potential energy; (c) (F c /F g ) eq , equato- 
rial ratio of centrifugal to gravitational force; (d) 
E, total energy, in units of GM 2 /(R e ) - The 
cooling evolution of a star that is represented by 
this sequence begins with the spherical model, for 
which K/Kq=1, and contracts along the sequence 
characterized by monotonically decreasing K/K . 
Quasistationary evolution proceeds until the bar 
mode instability is encountered at /3t ar ~ 0.27. 
The model for which jj=(j bar ~ 0.27 is marked 
with a solid dot in each frame. 




Fig. 2. — Density contours of selected models on 
the n'= 3 equilibrium sequence are shown in the 
(x > 0, z > 0) plane. The maximum density p max 
is normalized to unity. The highest density con- 
tour level is 0.9; subsequent contour levels range 
from 10 _1 to 10~ 10 and are separated by a decade. 
The axis ratios R p /R e of the models displayed are 
(a) 1.00; (b) 0.700; (c) 0.300; (d) 0.004; (e) 0.002. 
The model with i? p /i? e =0.004 shown in (d) has 
l3=/3 bar ~ 0.27. 
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Table 1 
n'=3 SEQUENCE. 



Rp/ Re 


K/K 




M 




J 




E 





(F c / Fg)eg 




VE 


1.000 


1 


,000 


7 


,73E- 


■2 


0.00 


0.264 


1.68E-8 


0.00 


0. 


,000 


1 


.87E-6 


0.900 





,989 


6 


.09E- 


■2 


5.06E-4 


0.223 


-3.16E-5 


5.28E-3 


0. 


199 


2 


.10E-6 


0.800 





,978 


4, 


,62E- 


■2 


4.61E-4 


0.183 


-4.01E-5 


1.07E-2 


0. 


,378 


2 


.58E-6 


0.700 





,967 


3 


.35E- 


•2 


3.37E-4 


0.146 


-3.53E-5 


1.62E-2 


0, 


,531 


3 


.46E-6 


0.600 


0. 


,955 


2 


.28E- 


■2 


2.11E-4 


0.112 


-2.52E-5 


2.19E-2 


0. 


,655 


3 


.94E-6 


0.500 


0. 


,942 


1 


.44E- 


•2 


1.13E-4 


8.11E-2 


-1.49E-5 


2.82E-2 


0. 


,749 


4 


.67E-6 


0.400 


0. 


,926 


8, 


13E- 


•3 


4.99E-5 


5.45E-2 


-7.13E-6 


3.54E-2 


0. 


,814 


6 


.67E-6 


0.300 


0. 


.906 


3 


.88E- 


■3 


1.69E-5 


3.26E-2 


-2.60E-6 


4.48E-2 


0. 


,854 


1 


.25E-5 


0.200 


0. 


,874 


1 


.39E- 


•3 


3.66E-6 


1.58E-2 


-6.04E-7 


5.91E-2 


0. 


,879 


2 


.56E-5 


0.100 





,807 


2 


,52E- 


•4 


2.86E-7 


4.69E-3 


-4.95E-8 


8.79E-2 


0. 


,897 


7. 


.84E-5 


0.050 


0. 


,725 


4 


.90E- 


•5 


2.47E-8 


1.41E-3 


-4.17E-9 


0.122 


0. 


,911 


2 


.28E-4 


0.040 


0. 


.696 


2 


.94E- 


•5 


1.15E-8 


9.64E-4 


-1.89E-9 


0.134 


0. 


,916 


3 


17E-4 


0.030 





,657 


1 


.53E- 


•5 


4.34E-9 


5.90E-4 


-6.86E-10 


0.150 


0. 


,922 


4 


.87E-4 


0.020 





,599 


6 


,24E- 


■6 


1.13E-9 


2.96E-4 


-1.66E-10 


0.174 


0. 


931 


8. 


.97E-4 


0.010 


0. 


,500 


1 


.41E- 


•6 


1.22E-10 


9.14E-5 


-1.49E-11 


0.215 


0, 


,945 


2 


.40E-3 


0.008 


0, 


,468 


8, 


.84E- 


•7 


6.09E-11 


6.28E-5 


-6.90E-12 


0.228 


0. 


,948 


3 


.18E-3 


0.006 


0. 


,430 


4 


.90E- 


•7 


2.52E-11 


3.89E-5 


-2.57E-12 


0.245 


0. 


,954 


5 


.29E-3 


0.005 


0. 


,406 


3 


.41E- 


•7 


1.47E-11 


2.88E-5 


-1.39E-12 


0.255 


0. 


,961 


7, 


.23E-3 


0.004 





,381 


2 


.21E- 


•7 


7.67E-12 


2.02E-5 


-6.62E-13 


0.267 


0, 


,964 


9 


.85E-3 


0.003 


0, 


,351 


1 


.31E- 


•7 


3.53E-12 


1.32E-5 


-2.68E-13 


0.280 


0. 


,973 


1 


.50E-2 


0.002 





,311 


7 


,20E- 


•8 


1.44E-12 


7.84E-6 


-9.66E-14 


0.297 


0. 


,992 


1 


.54E-2 
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